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Radio continuum surveys have, in the past, been of restricted use in cosmology. Most studies 
have concentrated on cross-correlations with the cosmic microwave background to detect the 
integrated Sachs-Wolfe effect, due to the large sky areas that can be surveyed. As we move into 
the SKA era, radio continuum surveys will have sufficient source density and sky area to play 
a major role in cosmology on the largest scales. In this chapter we summarise the experiments 
that can be carried out with the SKA as it is built up through the coming decade. We show that 
the SKA can play a unique role in constraining the non-Gaussianity parameter to ~ 1, and 
provide a unique handle on the systematics that inhibit weak lensing surveys. The SKA will also 
provide the necessary data to test the isotropy of the Universe at redshifts of order unity and thus 
evaluate the robustness of the cosmological principle. Thus, SKA continuum surveys will turn 
radio observations into a central probe of cosmological research in the coming decades. 
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1. Introduction 

Over the past decade it has become clear that large-area radio surveys can play an important role 
in enhancing our understanding of the cosmological model. Various groups have used the NRAO 
VLA Sky Survey (NVSS; Condon et al. 1998) as a low-redshift tracer of the large-scale structure 
of the Universe by both cross-correlations with the Cosmic Microwave Background (CMB) to find 
evidence for the Integrated Sachs-Wolfe effect (Sachs & Wolfe 1967; Giannantonio et al. 2008), 
and by investigating the largest scales for signs of non-Gaussianity (e.g. Xia et al. 2010). 

As we move into the SKA-era, then such wide-area surveys will reach much deeper levels, 
and begin to rival optical surveys in terms of source density. This is because the synchrotron 
emission that dominates the extragalactic radio background emission at < 10 GHz is emitted from 
all galaxies with any ongoing star formation or accretion activity. In addition, the wide area radio 
surveys will extend to substantially larger areas (~ 30,000 deg^) than many of the forthcoming 
optical surveys. 

Given that synchrotron radio emission from galaxies is unaffected by dust obscuration, and 
the k—corrections are relatively straightforward given the power-law nature of the synchrotron 
spectrum, then radio surveys offer a unique advantage for many cosmological applications. On 
the other hand, the generally featureless spectrum means that redshift information for individual 
sources can only be obtained with observations at other wavelengths, or through 21-cm spectral 
line observations (see e.g. Bull et al. 2015). 

Like optical surveys, radio surveys are made up of sources that trace the underlying density 
field very differenlly, i.e. from highly biased fracers such as powerful AGN fo low-bias populafions 
such as slar-forming galaxies. However, unlike fhe wide-held opfical surveys, radio confinuum 
emission is unaffecfed by dusf, and in fhe age of fhe SKA, fhe slar-forming galaxies, as well as fhe 
AGN, can be delected lo high redshifls. These sources can Ihen provide an imporlanf probe of fhe 
very largesl scales in fhe universe. 

In Ihis chapler we review fhe differenl experimenls and surveys lhal could be carried oul wilh 
fhe SKA in phase 1, and in fhe longer lerm wilh Ihe full SKA, for cosmology using radio conlinuum 
observalions. 

2. Assumptions 

We base our eslimales of Ihe redshift dislribulion of various source populafions on Ihe luminosity 
funclions lhal underpin Ihe semi-empirical exlragalaclic sky simulalions of Wilman el al. (2008, 
2010). These simulalions provide a very good descriplion of Ihe lalesl source counls from various 
deep held surveys wilh Ihe JVLA (e.g. Condon el al. 2012). Allhough modiftcalions may be re¬ 
quired lo accurately reproduce Ihe mosl recenl resulls from e.g. Herschel, Ihe general Irends and 
evolulion prescribed are relalively well malched lo our currenl underslanding, and Ihe exlrapola- 
lions lo flux-densily levels yel lo be reached in Ihe radio band are conslrained by observalions al a 
range of olher wavelenglhs. For example, in Ihe chapter on weak tensing al radio wavelenglhs and 
lesling Ihe foundalions of cosmology Brown el al. (2015) and Schwarz el al. (2015), bolh use an 
increase in Ihe source density of Ihe slar-forming populalion by a factor of 2.5, in order to malch Ihe 
lalesl source counls from deep fields. In Table 1 we presenl Ihe number counl for each radio source 
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Table 1: Expected number density (deg^^) of the various radio source populations in the simulations 
of Wilman et al. (2008) for a survey with detection threshold of Si ghz = 100 nJy {top-left), IflJy {top- 
right), 5/zJy {bottom-left) and 10/zJy {bottom-right). SFGs=star-forming galaxies; SBs=starburst galaxies; 
RQQs=radio-quiet quasars; FRIs=Fanaroff-Riley Class I sources; FRIIs=Fanaroff-Riley Class II sources. 
Note that the number density of SFGs and SBs in the table may be multiplied by a factor of 2.5 to match the 
latest observed source counts at the /zJy level (see Brown et al. 2015, for further details). 



100 nJy 

1 /ijy 1 

Redshift 



Wdeg^ 





V/deg^ 




SFGs 

SBs 

RQQs 

FRIs 

FRIIs 

Total 

SFGs 

SBs 

RQQs 

FRIs 

FRIIs 

Total 

0.0< z <0.5 

5892 

35 

4530 

358 

0.0 

10815 

3419 

37 

1710 

279 

0.0 

5445 

0.5<z<1.0 

15366 

191 

7093 

2833 

0.2 

25483 

7759 

260 

1964 

1004 

0.2 

10987 

1.0<z<1.5 

19353 

731 

4311 

3191 

0.8 

27587 

9181 

666 

1466 

857 

0.8 

12171 

1.5<z<2.0 

18102 

599 

2478 

1920 

1.2 

23100 

8042 

719 

967 

598 

1.2 

10327 

2.0<z<2.5 

16455 

1103 

1828 

1706 

1.0 

21093 

5844 

697 

664 

449 

1.0 

7655 

2.5<z<3.0 

13572 

888 

1104 

1198 

1.6 

16764 

4141 

525 

438 

328 

1.6 

5434 

3.0< z <3.5 

11582 

925 

787 

895 

0.8 

14190 

2970 

369 

281 

249 

0.8 

3870 

3.5<z<4.0 

10298 

887 

613 

848 

0.3 

12646 

2213 

302 

204 

209 

0.3 

2928 

4.0< z <4.5 

8144 

790 

414 

579 

0.3 

9927 

1676 

232 

145 

163 

0.3 

2216 

4.5<z<5.0 

7026 

649 

306 

481 

0.2 

8462 

1259 

184 

102 

130 

0.2 

1675 

5.0< z <5.5 

6037 

720 

233 

415 

0.0 

7405 

1005 

154 

75 

110 

0.0 

1344 

5.5< z <6.0 

5257 

638 

182 

319 

0.0 

6396 

803 

126 

54 

93 

0.0 

1076 

Total 

137084 

8156 

23879 

14743 

6.4 

183868 

48312 

4271 

8070 

4469 

6.4 

65128 


5 pJy 

10 ply 1 

Redshift 



Wdeg^ 





V/deg^ 




SFGs 

SBs 

RQQs 

FRIs 

FRIIs 

Total 

SFGs 

SBs 

RQQs 

FRIs 

FRIIs 

Total 

0.0< z <0.5 

1761 

40 

701 

140 

0.0 

2642 

1117 

40 

445 

99 

0.0 

1701 

0.5<z<1.0 

2894 

267 

913 

425 

0.2 

4499 

1594 

208 

636 

293 

0.2 

2731 

1.0<z<1.5 

2973 

368 

713 

351 

0.5 

4405 

1588 

225 

495 

239 

0.5 

2548 

1.5<z<2.0 

2350 

317 

468 

245 

1.2 

3381 

1206 

192 

309 

167 

1.2 

1875 

2.0<z<2.5 

1462 

218 

290 

184 

1.0 

2155 

693 

129 

179 

124 

1.0 

1126 

2.5<z<3.0 

930 

162 

174 

136 

1.2 

1403 

405 

95 

97 

92 

1.2 

690 

3.0< z <3.5 

605 

118 

103 

105 

0.7 

932 

239 

67 

53 

71 

0.7 

431 

3.5<z<4.0 

407 

91 

66 

84 

0.3 

648 

145 

50 

31 

56 

0.3 

282 

4.0< z <4.5 

274 

69 

39 

66 

0.4 

448 

87 

38 

17 

44 

0.4 

186 

4.5<z<5.0 

185 

53 

26 

53 

0.1 

317 

53 

28 

9 

35 

0.1 

125 

5.0< z <5.5 

130 

42 

17 

42 

0.0 

231 

34 

21 

5 

28 

0.0 

88 

5.5< z <6.0 

91 

35 

10 

36 

0.0 

172 

21 

17 

3 

24 

0.0 

65 

Total 

14062 

1780 

3520 

1867 

5.6 

21235 

7182 

1110 

2279 

1272 

5.6 

11849 


population for four flux-density thresholds at v = 1 GHz, whieh are representative of the typical 
surveys that could be conducted with the SKAl and SKA2. However, we note that moving the 
frequency to around 700 MHz would lead to an increase in the field-of-view by factor of ~ 2 and 
an increase in the source density that scales with the spectral index, i.e. Sy v“, with a = —0.7 
for a typical star-forming galaxy. 

3. Evolution of bias for radio sources 

One of the key measurements that can be extracted from surveys is knowledge of the bias of the 
population, i.e. how the light emitted traces the underlying dark-matter distribution. It will be 
very challenging to measure the evolution of the bias of radio sources from the large-scale surveys 
that will predominantly be used for measuring the cosmological signal, due to the lack of redshift 
information. However, assuming that the bias of a given source population is scale-invariant, then 
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Figure 1: Linear bias of near-infrared identified radio sources as a function of median redshift. Open circles 
correspond to the four independent redshift bins used while the filled circle is the bias for the full sample of 
766 radio sources with Ks < 23.5 and 4 > 90/zJy. Star symbols correspond to lower luminosity limits of 
10^^, 10^^'^ and 10^^ W Hz^' from low to high redshift. The dashed line shows the FR-11 bias adopted by 
(Wilman et al. 2008) in the SKADS simulations, and the diamond symbols show the expected bias based on 
the SKADS prescriptions (open) and with the FR-1 halo mass increased to lO'^ Mq (filled), matching the 
FR-IIs. Taken from Lindsay et al. (2014a). 


deeper surveys with a greater degree of multi-wavelength information can be used to measure both 
the redshift distribution and the bias of the radio sources, in conjunction with measuring the large- 
scale structure for the wide-area surveys. Such data will also be important for understanding the 
effects that baryons have on the underlying dark matter, and how this in turn affects the clustering 
signal. 

An initial attempt at constraining the bias of the radio source population was carried out by 
Lindsay et al. (2014a,b) who used a combination of wide-shallow and deep-narrow radio and 
optical/near-infrared survey data to determine the clustering of the radio source population in red¬ 
shift shells (Figure 1). They found that the dominant AGN population, namely the low-luminosity 
(FRI-type sources as defined in the SKADS simulations) appear to have a higher bias than cur¬ 
rently simulated. This provides an opportunity to use the multi-tracer technique to overcome the 
debilitating effects of cosmic variance on large-scale measurements of the power spectrum using 
various probes (Seljak 2009). 

However, such work is not limited to small scales, where the ancillary data provides the radial 
information. With the SKA we will also be in a position to cross-correlate the radio source popula¬ 
tion with CMB-lensing maps from current (e.g. Planck Collaboration et al. 2013a; Das et al. 2014) 
and future CMB facilities. Some work along these lines has already been performed (Sherwin et al. 
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2012; Geach et al. 2013), providing a unique measurement of the bias for high-redshift quasars. 

With information about the redshift distribution of the various sub-populations of radio sources 
(e.g. star-forming/AGN type) in any deep radio survey, then this information can inform the mea¬ 
surement of the bias from cross-correlating with CMB-lensing maps. This can provide a robust 
measurement of both the weighted bias of all populations, or if they can be morphologically distin¬ 
guished (see Makhatini et al. 2015), the bias and its evolution for each individual source population. 

4. Cosmological Probes 

4.1 The Angular Power Spectrum 

In the the absence of redshift information, the most straightforward experiment is to measure an¬ 
gular correlation function or power spectrum of the radio sources. Raccanelli et al. (2012) and 
Camera et al. (2012) have already shown that a measurement of the angular power spectrum, in 
combination with the CMB and supernovae la, can be a very useful tool to determine both the dark 
energy equation of state and departures from general relativity, even with the precursor surveys 
within reach of ASKAP, LOFAR and MeerKAT. However, to make radio surveys competitive with 
the largest surveys at other wavelengths then we require the combination of source density, sky 
area and morphological characterisation that is only feasible with the SKA. 

Recently, Ferramacho et al. (2014) demonstrated that even without redshift information for 
the individual sources, or a subset of sources, the wide-area radio continuum surveys can play a 
unique role in constraining the level of non-Gaussianity. Utilising the multi-tracer technique (Seljak 
2009), they showed that the different populations of radio sources, which trace the underlying dark- 
matter distribution with vastly different biases, can constrain the local non-Gaussian parameter fyi 
with uncertainty = 3.6 for a galaxy detection flux limit of 10 /iJy and = 2.2 for 1 /rJy 
(Figure 2). The former survey is within reach of SKAl, but requires good resolution in order to 
morphologically distinguish the different classes of radio source, i.e. FRI/FRII from star-forming 
galaxies and radio-quiet quasars. As shown in Makhatini et al. (2015) this sort of classification is 
possible to very high redshift with SKAl-MID. 

Therefore radio surveys with SKAl-MID, without any additional data, have the potential to 
constrain primordial non-Gaussianity to a factor of ~ 2 better than the present constraints obtained 
with Planck. This leads to the possibility of obtaining = 2.2 with SKA2 (Figure 2). 

Furthermore, the late time ISW effect allows us to strongly constrain a class of models that 
can have wq and Wa very similar to a quintessence model, or even to ACDM, but that allows for 
viscous pressure. The set of models that are called unified dark sector models or similar (e.g. a 
single viscous dark fluid) could in principle resemble LCDM in and wq, and produce a much 
larger ISW effect (see e.g. Li & Barrow 2009; Velten & Schwarz 2011). 

4.2 Cosmic Dipole, Isotropy and Homogeneity 

All-sky continuum surveys with the SKA will probe large angular scales at a median redshift 
of order unity. This will allow tests of fundamental assumptions of modern cosmology, especially 
of statistical isotropy. Combined with spectroscopic redshifts, either from SKA HI surveys, or by 
means of other instruments, statistical homogeneity of the cosmos will be tested. SKA will allow 
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0.001 0.01 0.1 1 10 

Detection limit (yU-Jy) 


Figure 2: Forecast constraints on obtained with the multi-tracer method as a function of the flux-density 
threshold used to detect galaxies. The various populations considered are FR-I and FR-II radio galaxies, 
radio-quiet quasars, star-forming galaxies and starburst galaxies, with different biases, as described in both 
Wilman et al. (2008) and Ferramacho et al. (2014). We present the results obtained using the full sample of 
objects with an averaged effective bias and those obtained using the combination of 3 populations of radio 
galaxies (where SRG, SB and RQQ correspond to one population group), using 4 populations (where only 
SFG and SB are undifferentiated) and with a selection of 5 populations for z < 1 and 4 populations for z > 1 
(again with undifferentiated SFG and SB). We also show the result for the ideal case where all 5 populations 
could be differentiated over the entire redshift range of the survey. The horizontal line represents the best 
constraint obtained by the Planck collaboration (Planck Collaboration et al. 2013c). Taken from Ferramacho 
et al. (2014). 


us to study scales at z ~ 1 that have not been in causal contact since the first horizon crossing during 
inflation, and therefore contain information that was frozen in during cosmological inflation. 

SKA all-sky surveys will allow the measurement of the cosmic radio dipole almost as precisely 
as the CMB dipole. SKAl will constrain the cosmic radio dipole direction with an accuracy better 
than 5 degrees (Fig. 3), and SKA2 within a degree (at 99 per cent C.L.). Compared to todays best 
estimate based on NVSS data, this will be an improvement of a factor of 100 in the accuracy of 
the cosmic radio dipole direction for SKAl. This measurement could firmly establish or refute the 
commonly adopted assumption that the CMB and the overall large-scale structure frames agree. 

The CMB exhibits unexpected features at the largest angular scales, among them a lack of 
angular correlation, alignments between the dipole, quadrupole and octopole, hemispherical asym¬ 
metry, a dipolar power modulation, and parity asymmetries (Planck Collaboration et al. 2013b; 
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Figure 3; Left; Accuracy (in per cent) of the measurement of the dipole amplitude as function of fractional 
error on flux density calibration on individual point sources (sigma). All points are based on 100 simulations. 
Right: Accuracy (in degrees) of the measurement of the dipole direction. The horizontal lines denote the 
error due to shot noise for a dipole estimate based on 8 x 10^ sources (SKA Early Science) and 2 x 10^ 
sources (SKAl). See Schwarz et al. (2015) for further details. 


Copi et al. 2013a,b). Understanding the statistical significance of these anomalies is crucial, as a 
lack of statistical isotropy or Gaussianity could rule out the standard cosmological model. 

The precision of these CMB measurements is limited by our understanding of the foregrounds, 
and observational uncertainties are already much smaller than the cosmic variance at such scales. 
Therefore, it is very difficult to identify the cause of these anomalies without an independent probe 
at the same scales. 

The angular two-point correlation at angles > 60 degrees from SKA continuum surveys, as 
well as the reconstruction and cross-correlation of low multipoles will offer further insight into 
these puzzles, see Schwarz et al. (2015). 

Supplemented with spectroscopic information, several tests of homogeneity can be envisaged. 
One possibility is to probe a consistency relation of Friedmann-Lemaitre models, which involves 
the dimensionless comoving distance, the Hubble rate, and their derivatives (Clarkson et al. 2008). 

Thus, studying the large-angular scales in SKA continuum surveys might help resolve the 
puzzle of CMB anomalies and test the cosmological principle. 

4.3 BAO and RSDs with spectroscopic follow-up of emission-line sources 

Although not a primary case for the SKA, radio continuum sources are, by definition, active and 
exhibit emission lines in their optical spectra. This property can be utilised to perform efficient 
follow-up spectroscopy with the next generation of optical and near-infrared spectrographs, in 
much the same way that the WiggleZ survey targeted UV-bright sources (e.g. Blake et al. 2012). 
Such follow-up would provide an efficient method for extracting precise redshift information for 
a very large number of sources over a large fraction of the sky. This could provide new and ro¬ 
bust measurements of the Baryon Acoustic Oscillations up to high-redshift (potentially targeting 
emission line objects at z > 2), as well as redshift-space distortions. 

Rather than a key driver for the SKA, this demonstrates a unique synergy, and something that 
is planned to be tested by combining the LOFAR extragalactic surveys with the new multi-object 
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spectrographs (MOS) that are currently in development (see Section 5). 

4.4 Integrated Sachs-Wolfe Effect 

The ISW effect (Sachs & Wolfe 1967) manifests itself in the correlation between large-scale struc¬ 
ture and CMB temperature. In an Einstein-de Sitter universe, the energy gained by a photon falling 
into a gravitational potential well is exactly cancelled out by the energy lost upon climbing out of 
the well. In a universe with a dark energy component or modified gravity to the same effect, the 
local gravitational potential varies with time, and potential wells are stretched throughout the cross¬ 
ing time of the photons. This disparity between potential on entry and exit imparts a net blueshift 
on the incident photon and, equivalently, an increase in photon temperature. Specifically, fhis effecf 
is imporfanf af lafe fimes (z<l), once fhe Universe had begun ifs fransifion towards a dominanf dark 
energy componenf and accelerafed expansion. The ISW signal is small, compared to fhe infrinsic 
femperafure anisofropies in fhe CMB, acfing on large scales where cosmic variance mosf affecfs 
CMB uncerfainfies, and so cross-correlation wifh local large-scale sfrucfure, wifh exfensive sky 
coverage, is required in order fo produce a significanf resulf. Giannanfonio ef al. (2012) discuss fhe 
current state of ISW measurements, and potential problems. While several detections of the ISW 
effect had been made previously, cross-correlating the CMB maps from the Wilkinson Microwave 
Anisotropy Probe (WMAP) with radio, infrared, optical, and X-ray surveys (for results using radio 
surveys see e.g. Nolta et al. 2004; Raccanelli et al. 2008), these were all at rather low significance. 
Giannanfonio ef al. (2008) reached an increased significance of ~ 4.5a by combining surveys fo 
develop a fuller cafalogue of local large-scale sfrucfure. Wifh fhe lafesf CMB measuremenfs from 
fhe Planck satellife, and corresponding all-sky surveys of massive, low-redshiff galaxies (and an 
undersfanding of fheir relafionship fo fhe underlying dark mafter disfribufion), fhe ISW could be a 
powerful cosmological fool (Figure 4; Raccanelli ef al. 2014). 

4.5 Cosmic Magnification Bias 

The presence of mass alfers fhe geodesic followed by lighf rays, causing fhe deflecfion of fhose 
rays along fhe line of sighf when passing by infervening large-scale sfrucfure. This causes distor¬ 
tions and magnification in fhe observed images of disfanf asfronomical sources, and is known as 
gravifafional lensing. The shape disforfion of fhe background galaxies (‘cosmic shear’) can be used 
fo consfrain some cosmological paramefers (see below, section 4.6), as can fhe magnificafion of 
background sources by foreground lensing strucfures (‘cosmic magnificafion’). 

This magnificafion is governed by fwo effecfs. The increased flux received from disfanf 
sources due fo lensing has fhe effecf of bringing info fhe survey sources which would ofherwise 
have been foo fainf fo be defecfed. The lensing also sfrefches fhe solid angle, reducing fhe apparenf 
surface densify of fhe lensed background sources. Gravifafional lensing should, fherefore, leave 
a cosmic magnificafion signal in fhe angular cross-correlafion funcfion of fwo samples of sources 
wifh non-overlapping redshiff disfribufions, since fhe foreground sample lenses fhe background 
sample. This signal, fhe magnificafion bias, is dependenf on fhe balance sfruck befween fhe loss of 
sources due fo dilufion, and fhe gain due fo flux magnificafion. 

Using fhe galaxy-quasar cross-correlafion funcfion, Scranton ef al. (2005) made an 8a de¬ 
fection of fhe cosmic magnificafion signal of quasars lensed by foreground galaxies, bofh selecfed 
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Figure 4: Constraints on (/nl) with integrated Sachs-Wolfe measurements using SKA continuum surveys. 
The effect of having several redshift bins for the measurements of the /nl parameter for different flux-density 
limits that will be possible for various stages of the SKA. For full details see Raccanelli et al. (2014). 

from the Sloan Digital Sky Survey (SDSS). They found that bright quasars, with steep source 
counts, exhibited an excess around foreground structure, and faint quasars, with shallow source 
counts were in deficit. Since this first success, Hildebrandt, van Waerbeke & Erben (2009) have 
made detections with normal galaxies from the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS), and Wang et al. (2011) have done likewise in the far-infrared with Herschel. 

Raccanelli et al. (2012) have shown that the cosmic magnification signal in SKA pathfind¬ 
ers can contribute to cosmological constraints; the cross-correlations envisioned there are between 
radio surveys (for the background sample) and optical surveys such as SDSS (for the foreground 
sample, and removal of radio foreground). With SKAl and SKA2, the continuum survey can pro¬ 
vide the background sample, while the HI survey can provide the foreground sample and removal 
of radio foreground (in combination with optical surveys such as Euclid and/or LSST). This will 
allow us to provide constraints on bias, dark energy and gravity parameters, with a method quite 
different from lensing autocorrelation of source shapes; this affords an important check on system¬ 
atic effects which could enter into shape measurement. 

4.6 Weak Lensing 

Weak gravitational shear is the coherent distortion in the shapes of distant galaxies due to the 
deflection of light rays by intervening mass distributions. Measurements of the effect on large 
scales is termed “cosmic shear” and has emerged as a powerful probe of late-time cosmology over 


■ Pathfinders 10 pjy 
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O SKA2.5pJy 
O SKA 1 pJy 

— Planck 
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3 low-z 
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the last 15 years (e.g. Heymans et al. 2012). Since gravitational leasing is sensitive to the total 
(i.e. dark plus baryonic) matter content of the Universe, it has great potential as a very robust 
cosmological probe, to a large degree insensitive to the complications of galaxy formation and 
galaxy bias. One of the most promising aspect of weak leasing measurements is their combination 
with redshift information: such measurements are then a sensitive probe of both the geometry of 
the Universe and the evolution of structure over the course of cosmic time. In turn, these latter 
effects are dependent on the nature of the dominant dark energy component in the Universe and/or 
on modifications to the theory of General Relativity on large scales. 

To date, the field of weak leasing has largely been fhe preserve of opfical surveys due fo fhe 
much larger number densifies of background galaxies achieved in such surveys. However, Ibis will 
change wifh fhe advenf of fhe Square Kilomefre Array, which will reach galaxy number densifies 
of up fo ~ 10 galaxies arcmin^^ wifh SKAl and ~ 75 galaxies arcmin^^ wifh SKA2 (Figure 5). 
In addifion, as described in fhe weak leasing chapfer in fhis volume (Brown el al. 2015), fhe radio 
offers fruly unique approaches fo measuring weak leasing fhaf are (i) nof available fo opfical surveys 
and (ii) pofenfially exfremely powerful in minimising fhe mosf worrying sysfemafic effecls in weak 
leasing cosmology. In particular, fhe PSF of an inlerferomeler is completely deferminisfic (af 
leasf al fhe mid frequencies of fhe SKA), and fhe inlrinsic synchrofron specfrum is smoolh and 
fealureless, fhus fhe colour-dependenf PSF should be far more slraighlforward fo accounf for, cf. 
broad-band opfical imaging. Moreover, using complemenlary radio shape informalion, along wifh 
polarisation informalion, which will be measured by defaulf in fhe SKA conlinuum surveys, will 
help fo confrol shape measuremenf sysfemafics and fhe impacl of inlrinsic alignmenls. For SKAl 
fhe radio continuum surveys will fherefore be very complemenlary fo fhe nexl generalion of opfical 
weak-lensing surveys from bofh fhe ground (e.g. DBS) and space (e.g. Euclid). 

In fhe longer lerm, surveys wifh SKA2 will extend fhe reach of weak tensing beyond lhal of 
fhe optical mega-surveys of LSST and Euclid. Firsl, fhe fypical redshifls probed by SKA surveys 
will go beyond lhal of opfical surveys. In addition, fhe size of survey which SKA will be capable 
of achieving will by ~ 30,000 deg^, Iwice lhal of Euclid and only rivalled (or complemenfed) by 
LSST (see Figure 5). 

However, as highlighled in (e.g. Palel el al. 2015) shape measuremenf al radio wavelenglhs 
need a similar effort fo lhal spenl in exlracfing shear eslimales from ground-based optical surveys. 

5. The multi-wavelength requirements 
5.1 Spectroscopic redshifts 

As wifh mosf cosmological measuremenfs, greafer accuracy can be achieved if higher precision 
redshifls can be obfained. In parlicular, using radio sources as largefs for nexl generalion MOS 
facilities could make such surveys highly efficienl due fo fhe higher likelihood of oblaining redshifls 
from emission lines. 

As defailed in Jarvis ef al. (2015) fhe fulure MOS facilities such as MOONS (Cirasuolo ef al. 
2012,2014), Prime Focus Speclrograph (PFS; Takada el al. 2014) and fhe Maunakea Specfroscopic 
Explore ^ (MSB) will provide a very good complemenl fo SKA surveys, albeil only over around 

'http://mse.cfht.hawaii.edu 
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Figure 5: Left panel: The redshift distribution of source galaxies for a Stz steradian weak lensing survey 
with the full SKA. Also shown is the redshift distribution for the 15,000 deg^ survey with Euclid . The n(z) 
extends to higher redshifts in the radio survey and probes a greater range of cosmic history. Right panel: The 
corresponding constraints on a 5-bin tomographic power spectrum analysis. For both experiments, an RMS 
dispersion in ellipticity measurements of /mis = 0.3 is assumed, and the tomographic bins have been chosen 
such that the bins are populated with equal numbers of galaxies. Open triangles denote 1(7 upper limits on 
a bandpower. Note that only the auto power spectra in each bin are displayed though much cosmological 
information will also be encoded in the cross-correlation spectra between the different z—bins. See Brown 
et al. (2015) for more information. 

half of the SKA sky for both PFS and MSB due to them being sited in Hawaii. For the wider 
surveys, 4MOST (de Jong et al. 2012) would provide a basis for obtaining redshifts for the brighter 
star-forming galaxies, predominantly in the low-redshift Universe. 

5.2 Photometric redshifts 

The bulk of the radio sources detected at these faint levels will be too faint at optical wavelengths 
to obtain spectroscopic redshifts. We will therefore be reliant on photometric redshifts. In the early 
phases this will be from surveys that are currently underway (e.g. McCracken et al. 2012; Jarvis 
et al. 2013; de Jong et al. 2013; Edge et al. 2013; Banerji et al. 2014) 

As we move to the full operation of SKAl then we should also have LSST (see Bacon et al. 
2015) and Euclid (see Kitching et al. 2015), which will provide very deep imaging from the g-band 
through to H-band across a large swathe of the southern sky. 

6. Survey requirements for SKAl 
6.1 Deep 

Although the main cosmological science will come from large-area surveys with the SKA, it is 
also important to conduct a deep survey with the best ancillary data available. As outlined in 
various chapters (e.g. Jarvis et al. 2015; McAlpine et al. 2015; Takahashi et al. 2015; Bacon et al. 
2015), deep fields covered with the LSST will be prime candidates for undertaking a deep radio¬ 
continuum survey. Such data will allow accurate photometric redshifts, and therefore the N{z) of 
the radio sources to be determined to a flux-density limit much fainter than possible over the wider 
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area. Using these data will act as a benchmark for the assumptions needed for fully exploiting the 
wide-area surveys. The key issue is that the cosmic volume surveyed should be representative and 
not suffer significantly from sample variance, which could act to skew the redshift distribution if 
large-scale structures dominate certain redshift slices. Moreover, such deep-fields provide us with 
the necessary information to understand how the radio sources trace the dark-matter distribution as 
we move from the single-halo term to the two-halo term at a few Mpcs. 

In order to reach the requisite depth, to gain insights into the source population as a function of 
redshift from morphological classification, and to conduct the most efficient deep survey, SKAl- 
MID would be the facility of choice for this tier. This is because SKAl-SUR has a field-of-view 
much larger than the individual deep-drilling fields fhaf will be fargefed by fhe LSST, which are 
of order 10 deg^. Furfhermore, fhe resolufion required, af fhe relafively low frequency fhaf such a 
survey needs fo be conducfed fo reach fhe requisife source densify/depfh, in order fo classify objecfs 
morphologically, would be exfremely difficulf fo achieve wifh SKAl-SUR (for defails see Jarvis 
el al. 2015). 

6.2 Medium-deep 

In practice, if will be necessary fo add fo fhe deep surveys fo overcome sample variance, parfic- 
ular for fhe mosf highly-biased populafions al moderate redshifls. However, fhe key driver for a 
medium-deep survey wifh SKAl is fo undertake fhe firsl wide-area weak lensing survey al radio 
wavelenglhs. As delailed in Brown el al. (2015), a survey covering 5000 deg^ fo 1 Jt Jy rms in 
Band 2, will provide excellenl constrainls on fhe power speclrum fhaf are very complementary to 
weak-lensing surveys being conducted at optical wavelengths on a similar timescale. 

6.3 Wide 

The key survey for measuring the largest scales in the Universe is obviously a survey that can cover 
as much cosmological volume as possible. First, given that distant radio sources can be found in 
relatively shallow surveys due to their extreme brightness, then area is initially more important than 
extreme depth, and conducting an all-sky survey very early in the expansion up to SKAl would be 
extremely beneficial for cosmology applications. 

In ferms of which elemenf of SKAl is fhe besf for such a survey, fhen fhe quesfion is more open 
fhan for fhe deep and medium-deep tiers, as fhe areal coverage benefifs from fhe large field-of-view 
offered by SKAl-SUR, allowing in principle fhe survey fo be conducted on a shorfer timescale. 
However, fhe resolufion of ~ 1 arcsec af 1 GHz for SKAl-SUR precludes fhe science fhaf requires 
fhe morphological characterisation of radio sources (see e.g. McAlpine el al. 2015; Makhafini el al. 
2015), which is needed fo exploif fhe mulli-lracer fechnique and overcome fhe inhibiting effecls of 
cosmic variance on fhe largesl scales (Section 4. 1), or fhe differenl redshifl disfribulions of differenl 
populafions. Given fhaf fhe cosmic variance is fhe dominant source of uncertainty on measurements 
of the power spectrum on the largest scales, then this is a crucial issue. 

Although the large scales can be probed efficiently by SKAl-SUR, the additional benefits of 
high-resolution imaging, which pave the way for using the multi-tracer method using SKAl-MID, 
provide a unique probe of the largest scales. Given that SKAl-MID is required to conduct the weak 
lensing survey, and the expansion to all sky for SKAl-MID does not present a large increase in 
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survey time eompared to SKAl-SUR (a SKAl-MID survey ean be earried out at lower frequeney, 
< 1 GHz, thus inereasing the size of the primary beam, while still retaining the requisite resolution 
for morphologieal eharaeterisation of the radio sourees), then it would make sense to eonduet even 
the widest of surveys with SKAl-MID. 

7. Towards SKA2 

As we move towards SKA2 then the ability to eonduet a world-leading weak lensing survey with 
eompletely different systematies to any optieal survey beeomes very eompelling. A souree density 
of ~ 75 sourees per aremin^ over 371 steradians is within reaeh (assuming a detention threshold 
of Si GHz > 100 nJy), eelipsing the souree density from eurrent surveys, sueh as the Dark Energy 
Survey (e.g. Banerji et al. 2014), and is higher density than Euclid will have. Furthermore, the PSF 
at radio wavelength is analytieally predietable, and the additional information from polarisation 
(and Hi rotation eurves) means that the major souree of uneertainty from intrinsie alignments may 
be overeome. Thus, SKA eontinuum surveys will turn radio observations into a eentral probe 
of eosmologieal researeh in the eoming deeades. As shown above, SKAl may already lead to 
transformational findings. 
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